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Abstract 

The X-ray sources observed in nearby galaxies show 
brightness distributions (log(iV(> S'))-log(5') curves) which 
can be described as single or broken power laws. Sin- 
gle power-law distributions are often found for sources 
in galaxies with vigorous ongoing star-formation activity, 
while broken power laws are found in elliptical galaxies 
and in bulges of spiral galaxies. The luminosity break can 
be caused by a population of X-ray binaries which con- 
tain a neutron star accreting at the Eddington limit or 
by aging of the X-ray binary population. We show that 
a simple birth-death model can reproduce single and bro- 
ken power-law log(iV(> S'))-log(S') curves. We have found 
that power-law log(A^(> S'))-log(5') curves are a conse- 
quence of smooth continuous formation of X-ray bina- 
ries and a luminosity break is the signature of a star- 
burst episode in the recent past. The luminosity break 
is robust and is determined by the lifespans of the X- 
ray binaries relative to the look-back time to the star- 
burst epoch. The model successfully explains the different 
forms of log(A^(> iS'))-log(5') curves for the disk and bulge 
sources in the spiral galaxy M81. 
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1. Introduction 

Recent Chandra observations have shown that nearby spi- 
ral galaxies similar to the Milky Way (e.g., MSl. lTennant et 
contain hundreds of discrete X-ray sources. The source 
numbers are even large in giant ellipticals, starburst galax- 
ies and mergers fsee lSarazin et al. 20001 IMatsumoto et al. 
IFabbiano et al. 200 While these sources are an inhomo- 
geneous population, the majority of them are believed to 
be X-ray binaries (XRBs) containing a black hole or neu- 
tron star. The X-rays from these binaries are powered by 
accretion, with the duration of their X-ray phases limited 
by the supply of accreting material from their compan- 
ion stars, which are mostly main-sequence or giant stars. 
The lifespans of XRBs are therefore finite. If the averaged 
mass-transfer rates are higher for brighter systems, then 
brighter XRBs are short-lived. 



For the X-ray luminosities of the Chandra sources in 
galaxies at distances of a few Mpc {L^ > 10'^^ erg s~^), the 
duration of their mass-transfer (X-ray) phase are shorter 
than 10^ yr. These timescales are short in comparison with 
the dynamical timescales of galaxies, which are usually 
longer than 10^ yr. XRBs are thus formed and then cease 
to be X-ray active within the lifetime of their host galax- 
ies. As the formation of XRBs is a ongoing process like the 
formation of normal stars, XRBs populate not just star- 
burst galaxies but galaxies of different ages and of different 
Hubble types. 

If a galaxy is undisturbed so that its stellar population 
ages smoothly, then the formation of XRBs is also likely 
to be a smooth process. If the galaxy has an episode of 
vigorous star-formation triggered by galactic interaction 
or merging, then enhancement in the rates of XRB forma- 
tion will subsequently occur. The XRB populations are 
therefore tracers of the dynamical history of galaxies. 

Here, we construct a birth-death model and use it to 
describe the time-dependent evolution of the XRB pop- 
ulations in galaxies. We show that the model is able to 
qualitatively explain the morphology of X-ray luminosity 
functions of XRBs in nearby galaxies. In particular, the 
model provides a mechanism for the formation of a lu- 
minosity break in the \og{N{> S'))-log(S') curves of the 
XRBs (where S is the observed X-ray flux). The implica- 
tions of the results from this study on our understanding 
of galactic dynamics and evolution are briefly discussed. 

2. Birth-death model 

For an XRB population in a galaxy with the number of 
within a luminosity range {L, L + AL) at time t 



specified by n{L,t)At, the cumulative luminosity distri- 
bution of XRBs is 



^>L,t)=y" dL'n{L', 



(1) 



As the sources are more or less at the same distance d, 
we can adopt A^(> S,t) = A^(> L,t), where the observed 
flux S oc L/And"^. For a population of systems determined 
only by birth and death processes (the beginning and end 
of the X-ray active phase, respectively), the evolution of 
n{L,t) is governed by a linear differential equation of the 
form 



f{L,t)-g{L,t) 



(2) 
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where f{L,t) is the birth rate and g{L,t) is the death rate. 

2.1. AN IDEALISED CASE 

As a simple illustration, consider the case in which (i) 
the sources have a finite X-ray phase duration ty^{L), say, 
proportional to 1/L, (ii) the source brightnesses are con- 
stant throughout the X-ray phase, and (iii) the sources are 
born at the same time, ta- We further consider a birth- 
rate function f{L,t) = \{L)5{t — t^), where \{L) is the 
brightness distribution of the sources at birth, and 5{..) 
is the usual Dirac (5-function. In an absolute deterministic 
case, the death-rate is a time translation of the birth rate: 
g{L, t) = /(L, t - t^{L)) = \{L)5{t -t^- t^L)). 
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Figure 1. The starhurst component of the N{> L) curves 
calculated from the idealised birth-death model (solid 
lines), without differentiating black-hole and neutron-star 
binaries. We have assum, ed that t^{L) = wh ere 

(3 = [ O.SMqC^ , and that that initial brightness distri- 
bution X{L) oc L~^-^, such that the N{> L) curves have 
slopes similar to that of the Chandra sources in M81 (see 
Tennant et al. 2001). Three cases, with {t — t^) = 0.1, 
0.5 and 1 Gyr, are shown. Their luminosity cutoffs are at 
1.7 X 10-''^, 3.4 X 10-''''' and 1.7 x lO'"'''' erg respectively. 
The N{> L) curve for the continuous star-formation com- 
ponent is also shown for comparison. We consider only the 
case in which the slope of hqIL) is the same as the slope 
of \{L) and the normalisation is 1/2 of the normalisation 
of the starburst component. 

When f{L,t) and g{L,t) are so specified, we can inte- 
grate equation (2) with respect to t, yielding 

n{L,t) = no(i) + XiL)9it - t^)[l - 9{t -t^- t^{L))] , (3) 

where 9{..) are the Heaviside unit step-functions. The con- 
stant fioiL) is zero, because of the assumption that no 
systems are born except at t^- The result here is clear — 
we only see XRB after a time t^, and we will not see 
systems with L after a time ta + t^{L). Moreover, the 
bright systems disappear first, followed by fainter systems. 



If t^{L) = {(3L)~^, where /3 is a constant, then a step-hke 
luminosity cutoff is formed at Lc = [P{t — ta)]~^ in n{L, t) 
and in the cumulative luminosity distribution iV(> L,t). 
As the XRB population ages, the cutoff migrates from 
high to low luminosities (Fig. 1). 

Galaxies also have continuous formation of XRBs due 
to ongoing star-formation (such as young high- mass XRBs 
in the arms of spiral galaxies) or orbital and nuclear evo- 
lution of binary systems (such as old low-mass XRBs in 
elliptical galaxies or in bulges of spiral galaxies) . This will 
produce a population of XRBs in addition to those formed 
as a result of starbursts. If t is large such that the birth 
and death of the XRBs due to continuous formation has 
reached an equilibrium, such that dn{L, t)/dt\con — 0, the 
linearity of the formulation allows us to simply set a non- 
zero no{L) in equation (3) and obtain a solution immedi- 
ately. Provided that the XRB population due to the star- 
burst dominates, the luminosity cutoff should be evident 
(Fig. 1). 

The formulation can be easily generalised to take into 
account the delay in the XRB formation after a starburst 
episode at t^. Wc can parametrise the delay e{L) as a 
function L, such that f{L,t) = X{L)S{t - (t^ + e(£))) 
and g{L,t) = \{L)5{t - (ta + e(i)) - t^{L)). It is again 
straightforward to integrate equations (2) and (1) directly 
and show qualitatively that the shape of the N{> L) dis- 
tribution of the XRBs obtained in the idealised case still 
hold. The only difference is that now the luminosity cutoff 
is no longer step-like. A more realistic and comprehensive 
modeling of the delay effects will be presented in Wu et 
al. (in preparation). 

2.2. A STOCHASTIC MODEL 

The purely deterministic approach given above needs to 
be modified because the onset of the X-ray phases of XRBs 
cannot be simultaneous and the systems may undergo 
transitions between bright and faint states. The observed 
luminosities L are not simply the lifespan-averaged lumi- 
nosities which can be related to the duration of the X-ray 
active phase. Also, the lifespans of the sources are charac- 
terised by their internal clocks, which are the individual 
source ages r, and they are not the same as the exter- 
nal clock by which the time parameter t specifying the 
evolution of the XRB population has been defined. The 
ages of individual systems are not direct observablcs in the 
study of populations of XRBs in galaxies. Only globally- 
averaged properties such as luminosity distributions can 
be measured. Hence, the death rate, defined in terms of t, 
cannot be considered as the time translation of the birth 
rate. 

If we assume that the systems with an observed lumi- 
nosity L have a death rate which is specified by a proba- 
bility mortality function k{L), and that the birth of XRBs 
is a Poisson process subject to a birth rate at time t, we 
can construct a simple statistical/stochastic birth-death 



formulation to describe the time-dependent evolution of 
the XRB population. The corresponding equation is 

j^n(L,t)^ f{L,t)-k{L)n{L,t) (4) 

IjWu 20011 IWu et al. 2002|l . As an approximation, we de- 
fine k{L) to be the inverse of the characteristic lifespan 
tx{L), such that t^{L) oc 1/L implies k{L) — PL, where 
P is a constant. For simplicity we neglect the delay in 
the formation of XRBs with respect to the starburst at 
time and assume that the birth function due to the 
starburst takes the form X{L)6(t — tg.). Thus, we have 
f{L,t) = Xo{L) + a\{L)5{t — t^), where Xo{L) is the con- 
tinuous component of the birth rate of XRBs and a is 
a parameter specifying the relative strength of the star- 
burst component. As f{L,t) and t^{L) are now specified, 
equations (4) and (1) can be solved, yielding n{L,t) and 
N{> L,t). 

An important property of the solutions is that, at any 
time t, the starburst component of n{L,t) has a cutoff 
at a critical luminosity Lc = [f3{t — ta)]~^. Here the cut- 
off is exponential, whereas in the idealised deterministic 
model considered above it is step-like. The location of the 
cutoff of the two cases is, however, identical. The cutoff 
causes an exponential luminosity break in the A'^(> L) 
distribution. Including time-delay effects in the formula- 
tion will not eliminate the luminosity cutoff but, instead, 
will cause A'^(> L) to appear as a broken power law for 
a power-law form X{L) (Soria, private communication). 
For the continuous component of n{L,t) there is no cut- 
off at Lc — [f3{t — ia)]~^- The corresponding cutoff is at 
Lc ^ [/3tgai]~^, where tgai is the age of the host galaxy. The 
luminosity break is therefore robust, and it is a character- 
istic of aging of a population of XRBs born as a conse- 
quence of an episode of star-formation in the recent past. 

We must emphasise that the results obtained from the 
above formulation should be interpreted in a probabilistic 
and statistical context. The formulation allows the pres- 
ence of systems with luminosity L brighter than the char- 
acteristic cutoff luminosity because the birth of XRBs is 
a Poisson process in the time domain and because the ob- 
served L is an instantaneous quantity while the duration 
of the X-ray phase of an XRBs is related only to the lumi- 
nosity averaged over the lifespan of the source. Moreover, 
n(L, t) is the expected value over an ensemble of sources 
with a distribution of ages which is not directly observ- 
able. (For more details, see IWu eF al. 2002 ) 

Nevertheless, the results obtained by the deterministic 
and stochastic formulations arc qualitatively similar, as 
can be demonstrated using the following example. Con- 
sider a simple birth-rate function, f{L,t) — X{L)[l+ad{t— 
^a)], where a is now the ratio of the strength of the impul- 
sive starburst component to that of the continuous forma- 
tion component. (Here we assume the same initial XRB 
brightness distributions for the two components.) Then we 



have 



n{L,t) = no{L) e'l^^' + 



ML) 

I3L 



1 



-l3Lt 



+aA(i)0(^-^a)e-^^^*-*^^ (5) 

where rioiL) = n(L,0) is the initial population. By com- 
paring the last term in equation (5) with equation (3) 
(with t^{L) = (fiL)~-^), we can easily see their similar 
time-dependent behaviours. Also, for PLt 1, the second 
term in equation (5) becomes X{L)/f3L, which is indepen- 
dent of t just as is the term tiq^L) in equation (3). 

Consider an initial population and a birth rate in power- 
law forms with the lower cutoff below the limit of detection 
i,, i.e., no = no,(L/L,)""i and X{L) = /o*(L/L*)""^ 
where ai and a2 are constants. Then the cumulative lu- 
minosity distribution is 

N(> L,t) = no,i,(/3L,i)"i-ir(l -ai,/?Li) 
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o-fo*Li, \fiL^, {t 



a2,(3Lt) 



X r{l~a2,PL{t-t,))e{t-t,) , (6) 

where r(Q;,a;) is the incomplete gamma function. 

Figure 2 shows four examples of iV(> L) calculated 
from this model. The luminosity break is clearly present 
for the impulsive-starburst component. By adjusting the 
relative strengths of the continuous and starburst compo- 
nents, we can obtain a form of A^(> L) which is a single 
power law or which has a luminosity break. Note the curve 
in panel (a) are very similar to what is observed in M81 
(Fig. 3). 

3. Discussion 

Above the completeness limit, the log(A^(> S'))-log(S') 
curves of XRBs in individual nearby galaxies can be de- 
scribed as a single or a broken power law. Separating the 
sources in the bulge and disk of the spiral galaxy M81 re- 
veals that the log(7V(> S'))~log(S') curves of the XRBs in 
the two galactic components are different IITennant et al. 2001|l . 
While the \og{N{> S'))-log(S') curve of the disk sources 
is a single power law, that of the bulge sources is a bro- 
ken power law. A power-law form for \og{N{> S'))-log(S') 
curve is generally found for XRBs in galaxies or galactic 
components with strong ongoing star-formation activity 
(e.g., the disk of MlOl. lPence et al. 20011 the nuclear re- 
gion of M83. Sori a & Wu 20021 and the starburst galaxy 
M82). Broken power-law log(A^(> S'))-log(S') curves are, 
on the other hand, found for the sources in early-type 
galaxies (e.g. NGC 4697,'Sarazin et al. 2000) and the bulges 
of spiral galaxies (e.g. M31, . Supper et al. 200 IJ 



Sarazin et al. (2000) proposed that the luminosity break 



is caused by a population of XRBs containing a neutron 
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Figure 2. The N(> L) curves of XRBs in a galaxy cal- 
culated from the stochastic birth- death model without dif- 
ferentiating black-hole and neutron-star binaries. The pa- 
rameters Ac = {fo*/no*PL^) and Xa = {afo*/no*) spec- 
ify the contributions of the continuous and the impulsive 
star-burst component respectively. The parameter f3 is set 
to be [0.3 MqC^]""'^, and the power-law indices ai and a2 
are both fixed at O.4. (0.): The starburst episode occurred 
4-00 million years ago (t ~ t^^ = 0.4 Gyr). The values for 
Ac and Aa are 1.0 and 0.23 respectively. The normalisa- 
tion is chosen such that the source counts are similar to 
the number of the M81 bulge sources observed by Chandra 
(see Tennant et al. 2001). The primordial component is 
represented by the dashed line; the continuous component 
by the dot- dashed line; and the impulsive starburst compo- 
nent by the dotted line. The total population is represented 
by the solid line, (b): Same as (a) except Aq — 0.03. (c): 
Same as (a) except Ac = 0.1. (d): Same as (a) except 
i - ta - 1 Gyr. 




star which accretes at the Eddington hmit. Using the 
stochastic formulation described above and assuming the 
Eddington luminosity limit (L ~ 2 x 10^^ erg s~^) for the 
neutron-star XRBs, it can be shown that a prominent lu- 
minosity break can be created fsee lWu et al. 20f)2)l . The 
break marks the Eddingtion limit for 1.5-Mq compact ob- 
jects. As this break is stationary and has the same value 
for all galaxies, it can be used as a distance indicator 
USarazin et al. 2000|l . 

While it is possible that the luminosity breaks in the 
log(A''(> S'))-log(S') curves of elliptical galaxies are caused 
by the luminosity limit of neutron-star XRBs, the lumi- 
nosity breaks in the log{A^(> S'))-log(5') curves of M31 
( |Shirey et al. 2001| ) and of the bulge sources in M81 IjTennant et al. 200l|l 
at about 4 x lO''^ erg s~^ must be due to another mecha- 
nism. 
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Figures. The (background substracted) log(N{> S))- 
log(S) curves of the disk and bulge sources in M81 ob- 
served by Ghandra. Adopted from Tennant et al. (2001). 



ulations of XRBs born at approximately the same time in 
the recent past. If this explanation is correct, then the 
location of the break marks the starburst epoch and the 
luminosity function becomes an indicator of the age of the 
host galaxy (or the host galactic component). The location 
of the break for the Chandra sources in the M81 bulge, for 
example, indicates that (i — ta) ~ 0.4 Gyr, which is consis- 
tent with the time M81 underwent a close encounter with 
its neighbor M82 ( |deGrijs et al. 2001| ) 
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From the birth-death model, we can see that a natural 
explanation for the luminosity break is the aging of a pop- 



